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For a lattice regularized chiral-invariant SU(2)l x SU(2)r fermion-scalar model with 
a Yukawa coupling y and a Wilson- Yukawa coupling w, we investigate the phase structure 
and in particular show the existence of the multicritical line, in the strong Yukawa and/or 
>- . Wilson Yukawa coupling region, at which four phases meet. The result is in good agreement 

/\ with the Monte Carlo simulation. This analytical result is derived from the effective scalar 

model obtained by integrating out the fermion field where the action is explicitly obtained 
from the hopping parameter expansion up to next-to- leading order. For estimates on the 
correlation function of the scalar field we apply the mean-field method. 
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In the standard model of electroweak interactions the Yukawa couplings of quarks and 
leptons to the Higgs scalar field provide the masses of the elementary fermions. In particu- 
lar, these masses are generated through the spontaneous symmetry breaking. However the 
formulation of the Yukawa models on the lattice suffers from the fermion doubling prob- 
lem [[|J]: the naively regularized Dirac fermion on the lattice describes additional fermion 
degrees of freedom (doubler fermions) besides one fermion desired in the continuum limit. 
To avoids this problem in the chiral theories one must remove unwanted doublers in the 
continuum limit in such a way that the chiral symmetry is preserved on the lattice [||, ||]. 
One formulation of the Yukawa model adds a manifestly chiral invariant Wilson- Yukawa 
coupling to decouple the doublers from the physical spectrum by rendering them heavy 
dynamically [|J]. The other method introduces mirror fermions explicitly, in addition to 
the original field, which transform oppositely under the chiral symmetry group [||]. 

In this letter we examine the former one, lattice Yukawa model with a chiral invariant 
Wilson- Yukawa term with coupling constant w, which reduces to the naive lattice Yukawa 
model with only a Yukawa coupling y in the limit w — > 0. The naive lattice Yukawa model 
(w = 0) which is defined below has a rich phase structure as shown in fig. [1], which is already 
discussed in the previous paper [||]. To distinguish the phase, two order paramters have 
been introduced: one is the usual expectation value of the scalar field (f> x and another is the 
expectation value of the staggered scalar field defined by (f> x := (— l) Xl + x ^ ^ Xd (f> x where d 
denotes the dimension of the Euclidean space. The phase diagram consists of the ferromag- 
netic (FM) phase (((f)) ^ 0, (0) = 0), the paramagnetic (PM) phase (((f)) = 0, (</>) = 0), 
the antiferromagnetic (AFM) phase (((f)) = 0, ((f)) 7^ 0) and the ferrimagnetic (FI) phase 
(((f)) 7^ 0, ((f)) 7^ 0). One of the most characteristic features of this phase diagram is that 
PM and AFM phases are separated into their weak (W) and strong coupling regions (S) 
by the presence of the intermediate FI phase. For w = the various phase transition lines 
intersect at the two quadruple points A and B where four phases can meet. 

At a nonzero value of the Wilson- Yukawa coupling 10/ the phase diagram can be 
approximately obtained by shifting in fig. [I] the y = axis by an amount of dw in the 
positive y-direction, as confirmed by various results of numerical simulations [0]. This is 
indeed understood from the leading result of the hopping parameter expansion given below. 
In the 3-dimensional bare parameter space (y,w,n) the two points A and B become two 
lines which we shall again label by A and B respectively. 

Quite recently Bock et al. [0] have shown that in the weak coupling region including 
the point A the unwanted fermion doublers can not be removed from the physical particle 
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spectrum. This implies that in the weak w region the continuum limit is not appropriate 
for describing the electroweak standard model. At least for w = 0, moreover, the existence 
of a nontrivial fixed point at A is quite unlikely from the result of the simulation for the 
renormalized Yukawa coupling [|J] . Thus the weak coupling region in the neighborhood 
of the origin y — w = is governed by the Gaussian fixed point up to the point A. For 
w — 0, Jiilich group [[JjJ] has also shown that the fermions become massive due to the strong 
Yukawa coupling y near the point B, which prevents the theory from having the meaningful 
continuum limit. Taking into account these facts, the physically meaningful point in the 
continuum limit is expected to reside in the vicinity of the multiple (would-be multicritical) 
point B in the weak naive- Yukawa coupling y <C 1 and strong Wilson- Yukawa coupling 
w ~ 1. 

This paper is the first to show by an analytical method that even in the presence 
of the Wilson- Yukawa coupling there are multicritical points, corresponding to B in the 
naive case w = 0, in the strong Yukawa and/or strong Wilson- Yukawa coupling region. 
This might corresponds to a nontrivial fixed point at which an interacting continuum limit 
can be realized. Moreover in the presence of the Wilson- Yukawa coupling the doubler mass 
can get mass of the order of the cutoff and therefore they are removed from the physical 
spectrum in the continuum limit. Basic ingredients are the hopping parameter expansion 
in higher orders and the estimates on the expectation value of the scalar fields, of which 
the mean field is the simplest one. 

The chiral SU(2)l x SU{2)r symmetric Yukawa model on the lattice is defined by 
the action S — Sh + Sp + Sy + Sw with 



Sh = ~kJ2 |tr($t$ x+M + $J$ x _^ 7 

X,fl 

X,fJ, 

x 

S w =w Y,{M®xPr + $IPl)1>x 

X,fJ, 

- ±M$xPr + *t +fl PL)ll>x +l * ~ &x($xPr + 4V^-M}> 

where k is the hopping parameter for the scalar field, y the normal Yukawa coupling, w the 
Wilson- Yukawa coupling, and Pl and Pr are the left- and right-handed chiral projectors, 



P L = |(l — 75), Pr = 5(1 + 75). Here the scalar field $ x is radially fixed $ x $ x = 1 
and is expressed by 577(2) matrix and the fermion field i\) and tfj are 577(2) doublets. The 
action has the global chiral SU(2)l X SU(2)r symmetry: ip x — >■ (QlPl + QrPr) il>x, 
4> x -» ^ (0[P L + OJjPfl), $ x -> L $ x O{?, where Lii? e SU(2) L)R respectively. 

It is convenient to transfer into the newly defined fermion field: i// = [y + 
d w ]-i/2 ($tp L + p fl ) ^ a , ) ?// = [ y + dw]' 1 / 2 $ x ($ X P R + P L ). Then the fermionic part of 
the action ([!]) is rewritten as follows: 

s F + s Y + s w = Y,W* -*£ ##«,»(*> ® ] W y , (2) 

x x,y 

where the matrix Q XjV is defined by 

d d 

Qx,y = Yl\- W ~ l^t^x+^L + P R )}Sy,x+n + ^ [w + ^(^ X ^P L + P R )]S V)X -^ (3) 

/x=l A*=l 

with a hopping parameter 

~ 2(y + dw)' ^ ' 

Here we notice that the fermionic part of the action has the bilinear form in the fermion 
field. After integrating out the fermion field, the model is converted into the effective 
scalar model with the action: 

S eff = S H -N f In DetM, (5) 

where Nf is the number of fermion species and the matrix M x>y is defined by 

M XtV ($,tf)=(l-KQ($,tf)) . (6) 

\ / x,y 

Now we perform the hopping parameter expansion (HPE) [|J]. The first step is to 
expand In Det M in powers of the hopping parameter K as follows: 

lnDetM x , y = Trln[l-KQ] x>y = -^— K 2n ^ tr[Q Xl>X2 ■ ■ -Q Xan , Xl ], (7) 

71=1 X 1 ,...,X2n 

where tr should be understood over the space of the Dirac spinor and the group 577(2). 
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Our result of the HPE up to 0(K 4 ) is 

InDet M = --K 2 D^hw 2 tr(l) - 2 tr{$i$ x+l 
- \k a D J2 (2w\7d - 4) tr(l) - 12w 2 tr[$t$ x+fJt ] ) 

X,fl 

- l^Dj^^l^+A^^] +2tr[$t$ x+M $t$ x+M ]) 

X,fl 

-\k 4 D J2 (8tr[$t$ x+ ^t +At+ ^ x+M ]_4tr[$t$ x+M $t +At+ ^ ;c+;y ] 



X,fI,V 

where D = 2 d ' 2 for even d. Here note that the trace is taken over the group G and that this 
result is independent of the choice of G. In the next step we specify the group: G = SU{2). 
The radially fixed scalar field $ x is a 2 x 2 SU(2) matrix and therefore can be written as 
$ = T a (/> a (a = 0, 1, 2, 3) and its conjugate as <&' = T a cf) a where 

T Q = (cr , icr m ), f a = (a , -icr m ), (a = 0, 1, 2, 3; m - 1, 2, 3) (9) 

with a unit 2x2 matrix &$ and Pauli matrices a m . By making use of the following 
identities for traces: 

\ t?[T a fP] = 6 a 'P, (10) 

l t Y [T a T f3 T p f a ] = 5 a ^S p ' a - S a ' p S^' a + 8 a ' a 8^' p , (11) 

final expression of InDet M is obtained as follows. 

InDet M = -DK 2 ]T[2«, 2 - 20^ +/ J 

X,fJ, 



- M 4 £(u; 4 (7d - 4) - 6«; 2 ^^ +M ) 



x,^ 



DK A J2 faX+^&L, - *%+»<%-» + <M+m<^+m - 2 J (12) 



x,/i 



- dk* J2 ( 4 ^c+,CC+^ - 2 ^ 



X,fl,V 



< Px4>x+ij, < Px+v ( Px+h+i> + 'Px'rx+p+v'Px+ii'Px+v < Px < Px+v < Px+p ( Px+ii+v 

Thus we have obtained the action of the effective scalar model up to 0(K A ). 
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In order to estimate the two- and the four-point correlation functions of the scalar 
field, we apply the simplest mean-field (MF) method. We introduce the external field 
H a (a = 0, 1, 2, 3) and define the expectation ( • ) H by 



)h 



<H 



'H 



n 



exp 



2> 



x,a 



n 



exp 



(E" 



(13) 



x I 



In this framework of MF theory, we can compute the correlation function explicitly, using 
the fact: 

(^) H = H a /4 + 0(H 3 ), (^ x ) H = 5^/4 + 0(H 2 ), (14) 

e.g., ($J$i$] +2 $i) = const. - ^- + 0(H 3 ). For the MF free energy given by 



F(H) = (S eS ) H + J2 H a m H - lnZ ff) 



(15) 



the stability-violating condition of the PM phase [0] -^p F(H)\ u=o = 0, yields the crit- 
ical surface separating the FM phase from the PM phase in the three-dimensional bare 
parameter space (y,w,K,) as follows. FM(S)-PM(S) critical surface: 



Kl 



d 



\-N f [ K 2 d + K 4 (6d 2 - 3d)w 2 + -K 4 (2d 2 - d) 



(16) 



where Nf = NfD = Nf2 d ' 2 . To obtain the critical surface separating the AFM phase 
from the PM phase, on the other hand, we have only to evaluate In Det M in terms of the 
staggered field $. Using the translational invariance, we obtain 



(InDet M) H = -]-K 2 DVd(w 2 tr(l) + /tr$J$i\ 
- -K^Dv{2d{7d - 4)w 4 tr(l) + 12d(2d - l)w 2 (tr i>J$i 



H 



f 8d(d- l)(trdj$i$l +2 $i\ -4d(d-l)/tr<ljl>ii>f +2 i> 2 



(17) 



J &>. <tJ 



sTA AT, 



+ 4d(tr$i ) $ l $[ +1 $ l ) +2d(tr& $ 1 ${ ) $ 1 ) ]+0(K 6 ), 

where V is the volume of the lattice. Applying the MF theory similarly, we obtain from 
the above condition AFM(S)-PM(S) critical surface: 



K2= d 



-1 - N f ( K 2 d + K*(Gd 2 - 3d)w 2 - -K i (2d 2 - d) 



(18) 



Now we pay attention to the intersection point of the two critical surfaces k = Ki, 
K = K2 just obtained. In the leading order (0(K 2 )) of HPE, there is no cross point as in 
the case of w = [||, flOU ]. In the next-to-leading order (0(K A )), the cross point does 



appear and makes a line in the bare parameter space (y, w, k), see fig. ^| (a). This is nothing 
but the multiple point at which four different phases meet. By equating both equation for 
the critical surface, k% — K2, the intersection line k = k c is obtained: 



y + dw 



N 



(2d 2 - d) 



1/4 



(19) 



which is depicted in fig. |2| (b). Indeed, in the limit w 

V4 



B = (yB,0,K B ): yB 



v,- 



(2d 2 - d) 



Kb — -'; 



2N f d 
2d-l 



0, this line has the end-point, 

-,1/2 

, which agrees with the 



previous result [||]. Another end-point in the limit y = locates at B' = (0, wg/ ,kb' 



w b > 



j_ 

2d 



N f 



(2d 2 - d) 



1/4 



K B ' 



2dN 



V,lV2 



2d-l 



{l + ^ (2d - 1)}, where w B > = 0. 



48, 



K B ' = -2.84 for N f = 4, d = 4. 

Although our explicit calculation in the presence of the Wilson- Yukawa term (w 7^ 0) 
is limited up to order if 4 , the existence of the multicritical line is guaranteed up to order 
0(K 8 ), since this line should have in the limit w = the cross point B whose existence is 
confirmed up to 0(y~ 8 ) in the previous paper [[f|]. Strictly speaking, we can not exclude 
the case that the multicritical line may terminate for large w and small y. However the 
Monte Carlo simulations [0] show that this is not the case. 

An important aspect of lattice construction which is not mentioned in this paper 
concerns the reflection positivity (RP). So far RP could be proven in a restricted range of 

@] 



the bare parameters in a restricted class of models [[ [IT 



. It is not yet clear whether 



RP holds also in the model treated in this paper [[0] and in particular in the region which 
we are focusing on. 

In conclusion, we have shown by the analytical method based on the hopping pa- 
rameter expansion up to next-to-leading order that the multicritical line at which four 
phases meet exist even in the presence of the Wilson- Yukawa coupling w for the chiral 
SU(2)l x SU(2)n symmetric lattice Yukawa model. This multicritical line is connected 
to the multicritical point B in the naive Yukawa model w = 0. In other words, the PM(W) 
phase and the PM(S) phase are separated by the intermediate FI phase. As a result it 
turns out that the PM(W) phase is not analytically connected to the PM(S) phase, and 
the physics in these two respective phases may be rather different. 
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In our model the scalar field is radially fixed which corresponds to the infinite self- 
coupling limit A — > oo of the scalar potential, A($$^ — l) 2 . It is a common wisdom that 
the radially fixed scalar model belongs to the same universality class as the weak-coupling 
one in the pure scalar sector. However recent numerical result based on large N shows 
that at small bare quartic self-couplling A the point A is replaced by a first order phase 
transition line, before the FM and AFM phases come close to each other [[[14[]. Therefore 



the multicritical point B' in the limit y — > is expected to be only one possible point for 
a nontrivial ultraviolet fixed point, since the vanishing of the fermion mass in lattice units 



in the limit is guaranteed from the Goiter man-Petcher symmetry [|15|1] and the doublers 
are removed by giving them masses of the order of the cutoff. 



There are however recent works [pj, [I7|, [0] which claim that the continuum limit 
K — > K c , y — > from the strong coupling regions FM(S) and PM(S) does most probably 
not lead to the desired formulation of the scalar-fermion sector of the Standard model 
on the lattice. It has been demonstrated in these works that the spectrum in the PM(S) 
phase consists only of the scalar particles ($) and a neutral fermion described by the field 
if;' . Moreover it has been argued that one ends up with an almost non-interacting neutral 
fermion in the continuum limit from the FM(S) phase Up. 61 1, since the renormalized Yukawa 



coupling of the neutral fermion to the scalar particles vanishes as a power of the lattice 
spacing a. In this scenario the dimensionless scalar expectation ($ ren ) goes to zero in the 
continuum limit. 

However we do not agree with this scenario and expect another one. In order to 
get the theory in the broken-symmetry phase by taking the scaling continuum limit, bare 
parameters must be adjusted in such a way that the vacuum expectation value of the 
renormalized scalar field ($ ren ) take a non-zero (finite) value in the limit. To perform this 
actually, we need detailed information on the phase structure and the critical behavior of 
the model as suggested from experiences in the scalar A^> 4 theory. The scalar Xcf> theory in 
the broken-symmetry phase constructed in the above sense has no self-interactions, i.e., the 
(dimensionless) renormalized three- and four-point self-couplings vanish in the continuum 
limit [[[19]]] • This is nothing but the triviality of the broken-symmetry \<p 4 theory. Thus 



we expect a scenario in which the non-vanishing Yukawa coupling is obtained by keeping 
the renormalized expectation value non-zero in the scaling limit by giving the triviality of 
the pure scalar sector as a sacrifice. We do not think that the triviality of the Yukawa 
coupling has been established at least in the case of approaching the line B-B' from the 



broken-symmetry phase. Detailed discussion on this issue will be given in a subsequent 
paper [ p0| . 



The behavior of the scalar propagator has already been investigated in the neighbor- 
hood of two multicritical points [[§§]. As a result, it has been shown that the scalar particle 
spectrum near the multicritical points A and B contains beside the usual bosonic particle 
additional 'staggered' particles which emerge as a consequence of the antiferromagnetic or- 
dering at negative k . The emergence of these spurious particles and the antiferromagnetic 
phases has to be simply regarded as an artifact of the special hypercubic lattice geometry 
and is a non-universal feature of this special kind of lattice regularization. However we 
do not think that the existence of two multicritical points A and B may be such a lattice 
artifact. This issue will deserve further studies. 
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Figure Captions 

Fig. 1. Phase diagram of the naive Yukawa model (w — 0) for Nf = 4, d = 4. 

Fig. 2. Phase diagram of the Yukawa model with the Wilson- Yukawa term for Nf = 4, 
d = 4: 

(a) schematic three-dimensional phase diagram, 

(b) projection of the phase diagram on the K=constant plane . 
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Fig. 2(a): Schematic three-dimensional phase diagram with the Wilson- Yukawa term. 
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Fig. 2(a): Schematic three-dimensional phase diagram with the Wilson- Yukawa term. 
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Fig. 2(b): Projection of the phase diagram on the K=constant plane. 
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